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A few subjects which strongly intertwine our field are discussed: K tttv decay, chiral symmetry on the lattice 
and a few other selected topics. Open questions are touched also on perturbation theory, locality, Gribov copies, 
CP symmetry in chiral gauge theories and cut-off effects. 



1. Introduction 

For the final talk of this exciting conference I 
have chosen a small number of subjects where the 
progress impressed me, or which dominated the 
field lately, or simply lie close to my heart. 

The first part deals with the heroic work of 
the CP-PACS and RBC groups on K ~* tttt. It 
is exciting and somewhat frightening to see how 
brains, software and hardware come together in 
a highly relevant project which is then pushed 
to the end (which does not necessarily mean full 
success). 

The second part is on chiral symmetric 
fermions. This field went through a rapid devel- 
opment during the last few years and attracted 
many contributions at this conference also. 

The third part contains selected topics: topo- 
logical susceptibility, some remarks on light 
hadron spectroscopy and a few words on cut-off 
effects. 

In a few cases I will raise questions, discuss 
issues which keep bugging me. These parts will 
be separated from the main text. It might be that 
some of them reflect my ignorance only. 



if ^ TTTT from K ^ t: and K 
domain v^rall fermions 
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There are many reasons for discussing this 
problem here. The K tttt decay contains highly 
relevant physics, it is a very active field of current 
research and a central problem of our commu- 
nity since the early years of lattice calculations . 
This complex problem is related to exciting theo- 
retical methods and issues like operator product 



expansion, (quenched) chiral perturbation theory 
{{Q)xPT), composite operator renormalization, 
chiral symmetry, etc. 

Interesting possibilities were discussed re- 
cently ||J^ to avoid the Maiani- Testa problem [Q 
due to three external particles in the K — > tttt 
amplitude. Preliminary results using three-point 
functions have been presented in Martinelli's ple- 
nary talk|^,^. Here I will discuss the 'reduc- 
tion method' used by the CP-PACS § and RBC 
groups ||]. In this case the K ^ mr matrixele- 
ment is related to the K n and K matrix- 
elements Without a chiral symmetric regular- 
ization the problem is not tractable. Both groups 
used domain wall fermions pO[|. 

The K — > TTTT decay exhibits two significant 
phenomena of the standard model: the A/ —1/2 
rule and the mixing induced and direct CP vio- 
lation. The dynamics of these decays is deter- 
mined by a non-trivial interplay of strong and 
electroweak forces characterized by several energy 
scales of very different magnitude from mt,mw 
to rrid, rriu- At present there is no way to treat all 
these scales at the same time on the lattice. To 
disentangle long and short distance contributions 
operator product expansion is used. 

2.1. Operator product expansion (OPE) 



Consider the 



TT transition in leading 



order of the weak interactions. The quark level 
transition on the l.h.s. of Fig. |l] is dressed with 
QCD interactions (the spectator quark is sup- 
pressed). Since Mw is large relative to the typical 
momentum p of this problem and so the W boson 
propagates over a short distance, the four fermion 
interaction is almost point like. This makes a sys- 
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(3-loop) PT up to as F ~ 0.3 in pure Yang-Mills the- 
ory [nd|. Having the A-parameter from this work (in 
terms of the low energy scale ro « 0.5fm[Q), PT in 
Yang-Mills theory is a parameter free expansion. In 



Figure 1. The leading order quark level graph 
and the corresponding four fermion interaction. 



tematic expansion possible for this amplitude: 

A = C{^i/Mw,as){Q) (1) 

with 0{p'^ /AI^) corrections. QCD effects above 
the factorization scale ^ are included in the Wil- 
son coefficient C, the low energy contributions be- 
low n are collected into the matrixelement of the 
local four fermion operators Q. This concept of 
factorization lies at the hart of many QCD appli- 
cations and I am not aware of reasons to question 
its applicability. 

If ^ Aqcd, we expect that the Wilson co- 
efficient C(/i/Mvi/, as) can be calculated in per- 
turbation theory (PT). The factorization implies 
that C is independent of the external states and 
can be calculated using off-shell quarks with high 
virtuality on both sides of the equation in Fig. |^. 
The scale fi has to be chosen judiciously. It should 
be large enough to justify PT when calculating C. 
On the other hand, if only three light quarks are 
simulated on the lattice (charm is integrated out), 
fi can not be larger than rric ~ 1.3 GeV. 

While factorization is a basic assumed feature, the 
question whether the Wilson coefficients can be calcu- 
lated perturbatively at a given n/ A.qcd is a technical, 
but important issue. 

One can check consistency within PT itself by con- 
sidering the convergence properties of the expansion, 
but this is not easy and also not the whole story. 
With the lattice formulation we have a tool to check 
the performance of PT by comparing it with the full 
non-perturbatively calculated result. PT is used in 
different situations: 

i) Physical quantities with a high scale 

The number of available well measured short distance 
physical quantities is very limited. The scale depen- 
dence of the running coupling in the Schrodinger func- 
tional scheme asF(M) is in good agreement with 
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Figure 2. The measured static force compared 
with PT (solid line) and the string model predic- 
tion (dotted line)0. 



a recent work[|lJ| the potential/force was measured at 
short distances 0.05fm < r < O.Sfm deep in the con- 
tinuum limit and the result was compared with PT. 
The authors pointed out that the perturbative scheme 
adopted should be chosen judiciously - a point which 
sheds additional light on the troubled history of PT 
for the quark potential^]. The scheme defined by 
the force aqq{fi) = F{r)/CF, fJ- = i/r looks well be- 
having and describes reasonably the non-perturbative 
data up to r x 0.15fm « (1.3GeV)"\ Fig. | 
shows another well known feature which is, how- 
ever, difficult to understand. Including the leading 
bosonic string model correction in the force F{r) = 
a + 7r/12r^||l^, which is expected to work at large 
distances, describes the data well at quite short dis- 
tances also. Actually, even at r ~ 0.15fm this form 
is closer to the data than PT. The picture is confus- 
ing since, even at the largest distance of this analysis 
(1.6fm) the spectrum of the effective QCD string de- 
viates significantly from that of an effective bosonic 
stringjl7[. 

ii) Renormalization constants Z 

The scale dependent renormalization constants of 
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Figure 3. The Z- factor of the vector current | 



composite operators depend on the cut-off and on 
the renormalization point /j.. When combined with 
renormalization group (RG) this should be a con- 
trolled perturbative problem if these scales are large. 
Nevertheless, as it is well known, PT performs poorly 
on the renormalization factors. The non-perturbative 
{1-loop boosted perturbative} pseudoscalar renor- 
malization factor Zp, for example, at fi = 2GeV, 
1/a « 2GeV is 0.45(6){0.62}@, 0.39(3) {0.59} 
and 0.34{0.54}[^ for Wilson, non-perturbatively 
improved Wilson and staggered actions, respec- 
tively [^. The situation is not always better with 
the finite renormalization factors, like that of the 
naive vector current Zv- This renormalization factor 
has no scale {fi) dependence, it is a function of the 
lattice coupling go and goes to 1 in the continuum 
limit. Fig. taken fromjiij ^, contains quenched 
a G 0.11 - 0.20fm and iV/ = 2 a £ 0.11 - 0.23fm data 
(RG improved gauge action , mean field improved 
clover quark action) compared with bare and mean- 
field improved PT. 

iii) Cut-off effects, improvement coefficients 

In order to eliminate the leading cut-off effects (in 
spectral quantities) the action is extended by a term 
Csw{go)Q, where Q is the Sheikholeslami-Wohlert 
operator j2^, while Csw is the Symanzik improve- 
ment coefficiently. The improvement should work 
for any external states, also for those with high vir- 
tuality, which can be used to calculate the improve- 
ment coefficients Csw (go)- (Notice the analogies 
with OPE.) For cut-off values of present day simula- 
tions (1/a < 3GeV) PT does not reproduce the non- 



perturbatively fixed (751^(50)1^71. For further exam- 
ples of problems with PT when calculating operator 
improvement coefficients see Table I in|28|. 

The conclusion is that in the very few cases where 
PT was systematically compared with numerical re- 
sults on physical quantities a judiciously chosen PT 
scheme seems to work well down to such low scales as 
fi ~ 1.3GeV. On the other hand, in problems where 
the cut-off also plays a role PT performs poorly for 
1/a < SGeV. Many believe that this is related to 
some special feature of lattice PT, like the existence 
of tadpole graphs. I do not think that this is the 
full explanation. First, a lot of effort was invested 
to overcome this problem of lattice PT[^. Second, 
more importantly, some of the quantities in question, 
do not look like perturbative at all. The a dependence 
of the coupling go of the plaquette action (lattice /3- 
function)[Q, the go dependence of Csvi/|2^], or Zv 
in Fig. 1^, for example, simply do not invite to use PT 
at all. 

The graph in Fig. |^ actually generates two four- 
fermion operators with different color structure 



Ql = {SaUb)L{ubda)L , 
Q2 = [SaUa)L{ubdb)L , 



(2) 



where a, b are color indices and the short-hand 
notation {saUb)R/L = Sa 7^(1 =•= 75)^6 is used here 
and below. The sum of these operators (weighted 
with the Wilson coefficients which include Vus KTd 
from the CKM matrix) taken between the incom- 
ing kaon and the outgoing pions in isospin / = 0, 
or 2 states 



K° (tttt)/ = Ai exp{i5i) = 

2 

should explain the large ratio 
1 Re(Ao) 



(3) 



Re(A2) 



22. 



(4) 



am indebted to Tomotcru YosMe for ; 
concerning the results on Zy in[|22| andl 



correspondence 



i.e. the A = 1/2 rule. In eq (y) Sj is the phase 
shift from the ffiial state tttt interaction, and 
Zi(/i) c>c Ci{ii) Wilson coefficients (and the addi- 
tional ones for i = 3, . . . , 10 below) are known in 
the standard model at next-to- leading order |pi[. 
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On the other hand, the weak interaction part of 
the 'tree graph' in Fig. |l| does not know about 
the third generation and is so CP conserving. In 
order to see CP violation we have to go beyond 
this leading approximation and include the QCD 
penguin and electroweak diagrams of Fig. ^. 
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Figure 4. The QCD and electroweak penguins. 



The resulting effective hamiltonian is a sum over 
10 operators [in. The A = 1/2 proble m IS ex- 
pected to be dominated by the Qi and Q2 contri- 
butions since the fimctions Zi{fi = rric) are small 
for i = 3, . . . , 10 ||3l[]. On the other hand, phe- 
nomenological considerations suggest [Q that the 
imaginary part of the amplitude 



ImAj = 



(5) 



V l^) .,=3 

where (Qi)/ exp(zi5/) = {Tnrj\Qi\K), is dominated 
by (Q6)/=o and {Qs)i=2- In this approximation 
the direct CP violation is given by 

e' _ ljGf 

7 ~ 2|e|ReAo' 



10 



-Im(T4,y,*) X 



(6) 



UJ 



where 

Qe = {sd)L {iuu)R + {dd)R + (ss)_r) , (7) 
Qs = ^{sd)L {('2uu)r - {dd)]i - {ss)r) . 

It is reassuring that this qualitative expectation 
is supported by the numerical data|0,^. Al- 
though the contribution from Qs is suppressed by 
the electroweak coupling constant, this is largely 
compensated by the prefactor 1/uj leading to a 
numerically unpleasant cancellation between the 
two terms in eq. (o). 



2.2. Determination of the bare 

{Tnri\Qi\K^) matrixelements 

The CP-PACS and RBC collaborations fol- 
lowed the 'reduction method', where tree level 
xPT is used to relate the K'^ tttt matrixele- 
ments to those of K'^ — > 7r+ and K vacuum^ 
which are then calculated on the lattice. 

Although the reduction step simplifies the 
problem significantly, it introduces an unpleasant 
mixing with a lower dimensional two-quark oper- 
ator. The {Sl,1r) operators, like Qe in eq. (0), 
mix with the two-quark operator 'sd' since the 
singlet right handed part communicates with the 
vacuum Q. This (8l, 1r) operator should be CPS 
(CP + d ^ s) invariant (as all the Qi operators 
are) leading to the form for md 

Qsub = {rus + md)sd - (tos - md)s^^d (8) 

'ms + md_ , ms-md_ , 
Oa I s7/^a ; s7m75« 



rus - md 



Wis + md 



In eq. the equation of motion was used to 
write Qsub as a total divergence for TOs 7^ m^. 
In the physical K'^ — > tttt matrixelement this to- 
tal divergence does not contribute since the weak 
operator carries zero momentum. On the other 
hand, for the —> tt"*" reduced matrixelement 
Qsub gives a (quadratically divergent) contribu- 
tion: if ms = md (jtItt = m^) is used in the sim- 
ulation (CP-PACS, 0) then Qsub is not a total 
divergence; if TOs 7^ md (RBC, [||) then the oper- 
ator carries momentum and the total divergence 
contributes. This quadratic divergence should be 
subtracted from the A/ — 1/2 operator matrix- 
elements with a coefficient determined by the cor- 
responding K —f matrixelement. The reduction 
equations have the form (/^ « 93 MeV) 



{Tr+TT-\Q 



A/=l/2 



(9) 



aiQsub\K^ 



2 9 

for i = 1, . . . , 6, 9, 10 and is determined by 

{0\Qf''='^^-a,Qsub\K+)=0. (10) 

The A/ = 3/2 operators have no subtractions 
and for z = 1, . . . , 6, 9, 10 the relations arc like in 
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eq. 



§) with 



0. For i = 7, 8 we have 



n+n-\Qi\K') = - 



V2U 



{n+\Ql\K^ 



In the equations above rn/f and TOtt are the phys- 
ical masses, while in the K —t n matrixclcnient 
TTiK = w,r = M . Remarks: 

• The reduction is on the tree level of (Q)xPT 
only. The 0(/) low energy constants needed for 
K —t TTTT can not be found from the K tt, 
K ^ Q matrixelements. One can use the O(p^) 
(Q)xPT results for if ^ tt |Q to make a more 
reliable fit for the tree level constants, but even 
in this case K —f tttt will be obtained on the tree 
level only. The 1-loop chiral corrections are esti- 
mated to be significant (m^ « 500MeV !)||,||. 

• The reduction step relies heavily on chiral sym- 
metry. The mixing problem can not be treated in 
the A/ —1/2 sector with Wilson fermions. 

• Both groups demonstrated that a controlled sig- 
nal can be obtained after the power divergent sub- 
traction. 

• As eq. (H) shows the if + 7r+ matrixele- 
ments should vanish in the chiral limit Af — s- 
for i ^7,8. CP-PACS finds that this is the case 
within the errors, RBC observes small non-zero 
intercepts which might be related to the resid- 
ual chiral symmetry breaking of the domain wall 
fermions. As eq. (H) shows, the slope is all what 
one needs. 

2.3. Renormalization 

CP-PACS used renormalization factors cal- 
culated in 1-loop PT Hi. RBC apphed the 
Roma-Southampton non-perturbative method 
with gauge fixing] 35 1. 



The uncontrolled effect of Gribov copies in the 
gauge fixing is a frequently discussed issue in this non- 
perturbative renormalization method. Undoubtedly, 
it would be very difficult to perform an unbiased aver- 
aging over the Gribov copies[^. On the other hand, 
is such averaging really needed if we study a short 
distance (oc l/fi) problem? Is it consistent to worry 
about Gribov copies and, at the same time accept- 
ing the validity of PT at the scale fi (say, in MS), 
which knows nothing about such non-perturbative 
features? 



2.4. Physical results 

The parameters of the simulations of CP-PACS 
and RBC are quite similar. On the real part of 
Aj and so on the A/ =1/2 rule, RBC obtained 
results close to those in experiments, while the 
results on KeAo and uj^^ from CP-PACS are a 
factor of ~ 2 smaller. 

Both groups obtained a small negative 0(10^"') 
number for e'/e as opposed to the experimental 
average of (17.2 ± 1.8)10-4[|7). 

Given the large theoretical and numerical com- 
plexity of the problem I find these works very im- 
pressive. There are several possible sources of sys- 
tematical error, like quenching, large corrections 
to tree level reduction in xPT and (continuum) 
PT at scales ~ l.SGeV. In[Q the authors argue 
that a new low energy constant (a pure quench- 
ing artifact) enters the reduction relations which 
might have an infiuence on the results above. 

3. Chiral symmetric lattice fermions 

We had four plenary talks at this conference 
which were directly on this subject js^, or were 
closely related ML I would like to discuss here a 
few points only|ll| . 

The Ginsparg- Wilson (GW) relation|Q was 
suggested as the 'mildest way' to break chiral 
symmetry on the lattice in 1982 immediately af- 
ter the no-go theorem of Nielsen and Ninomiya 
Bs]. In its simplest form the GW relation reads 



(11) 



or, equivalently (if D has no zero modes) 

-f5D~^{x,x') + D^'^{x,x')j5 = 754,x' • (12) 

Although the propagator D^^ does not anticom- 
mute with 75 the violation is a contact ~ S^^x' 
term only. It is expected and is really so that the 
GW relation implies chiral symmetry on physical 
predictions Q]. 

Eq. ( |ll| ) is a non- linear relation for D. It 
is clear intuitively and has been shown rigor- 
ouslyl^^ that the solution can not be ultralocal, 
it has a tail. On the other hand, a physically in- 
teresting solution should be local, which means 
that this tail of D should decay faster than the 
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signal at physical distances in correlation func- 
tions 

D{x, x') - exp {-Ma\x - x'\)^^_.^,^^^ , (13) 
where M = O (cut-off). 

Let me discuss the important issue of locality of 
the action in a QFT further. Locality is related to 
universality and universality implies predictive power. 
If, however, the action has couplings between distant 
points which compete with the real dynamics of the 
system, then the predictions will depend on the non- 
physical microscopical details. Universality will be 
lost and such an action is not useful. 

For this reason, the checks concerning the local- 
ity of the GW type Dirac operators are very im- 
portant | |39| , ^6[ -|49| . Dynamical staggered fermions 
with Nf < 4 present a not fully understood, po- 
tentially dangerous situation. It is not obvious that 
{detD(U)st)^^^ in the path integral defines a local 
theory, although there are no arguments excluding 
this either |5^. Two dimensional models might be the 
right place for testing. 



Since no solution to the GW relation eq. ( |ll| ) 
was found in the presence of gauge fields, the idea 
was soon after its birth abandoned. The first lat- 
tice regularization of fermions with chiral symme- 
try, the domain wall fermions and the related 
overlap constructionist | followed a different path 
and seemed to be unrelated to the GW relation. 

Following the observation that the fixed-point 
Dirac operator satisfies the GW relation|52| the 
interest turned to this general formulation again. 
The GW relation is a powerful theoretical tool. 
Its immediate consequence is the index theorem 
on the lattice |5^] and it implies the existence 
of an exact chiral symmetry transformation [Q. 
The observation that the overlap Dirac operator 
satisfies the GW relation connected the domain 
wall approach to this general formulation [55|] . 

As the previously discussed K — > tttt studies 
illustrate, domain wall and, to a lesser extent, 
overlap fermions are used today in large scale 
quenched simulations in problems, where chiral 
symmetry plays an important rolejs^J^,^. 

The fixed-point gauge and Dirac actions are de- 
fined by classical equations |p7[ . These equations 



can be solved numerically in an iterative way 
which is, however, far too expensive in a stochas- 
tic calculation. A parametrization in terms of a 
finite number of operators is unavoidable. This 
can be done quite effectively and the Dirac op- 
erator obtained this way performs well in differ- 
ent test runs ||5|,Q. An alternative way to fix 
the couplings in a parametrized Dirac operator is 
a systematic expansion of the GW equation it- 
self!^. For test results on this action and for 
comparison with other formulations, see[|60|. 

The new chiral symmetric actions allow to in- 
vestigate the role of topological excitations in the 
QCD vacuum. It was convincingly demonstrated 
in several contributions at this conference that an 
instanton dominated picture of low eigenmodes 
of the Dirac operator is consistent: the local 
peaks in these eigenmodes are dominantly chi- 
ral 1 61 6^. Such a study is hardly possible with 
Wilson fermions [Q. 

Chiral symmetric Dirac operators have nice 
theoretical properties including the absence of ex- 
ceptional configurations and 0{a) improvement 
in spectroscopy. They are, however, expensive to 
simulate. The first studies on spectroscopy are 
exploratory ||5^ , |63| - |65| , |58|| . From these works it is 
obvious that light hadron spectroscopy becomes 
much cleaner: quark mass tuning and long chiral 
extrapolations are not needed, topology is well 
defined and operator renormalization is signifi- 
cantly simplified. The first quantitative results 
confirm that there are no unexpected, hidden 
problems with this regularization of QCD. Simu- 
lations in the Schwinger model strengthen this 
conclusion. 

There are problems, like the fermion conden- 
sate, where a chiral symmetric action is a prereq- 
uisite for a quantitative study. The quark con- 
densate, or more precisely, the leading low energy 
constant S in the chiral Lagrangian, has been cal- 
culated by several groups using different actions 
and methods |63|6|,||1. The results are consistent 
and scatter around (I]^(2GeV))i/3 = 270MeV 
with a few percent statistical error. Since no con- 
tinuum extrapolation is done yet, the systemati- 
cal error is unknown. 

As mentioned before, eq. (O) implies exact 
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symmetry under the transformation |54[| 

^ iP + ieT^jsil - D)il; , (14) 

a = 1 , . . . , N'j — 1 , which is a sort of smeared 
75 transformation. The action is also invariant 
under a singlet chiral transformation, but the 
fermion measure is not (notice that the trans- 
formation in eq. ( p^ ) is gauge field dependent) 
leading to the correct anomalous Ward identity 

{60)F~2Nfv{U){0)F=Q. (15) 

where v{U) is the topological charge, while { ■) f 
is the fermionic path integral on a fixed gauge 
field configuration U . 

The transformation in eq. ( p^ ) is asymmetric 
on ip and "0, and so are the left/right projectors 
also Q. 

V'L/fl = , 4>L/R = i^PR/L , (16) 

where 

Pfl/i = i(l±75), Ph/l = ^(1±75), 

75 = 75(1 -i?). (17) 

The fermion action falls into L and R parts as in 
the formal continuum 

0-D0 = ^lD^l + ^rD^jr . (18) 

Eqs. ( p^ , p^ , p^ are the first, almost trivial steps 
in a highly complex construction which led to 
a breakthrough concerning chiral gauge theo- 
ries For anomaly free complex representa- 
tions for U{1), or SU{2) x U{1) gauge groups the 
theory is constructed non-perturbatively ||69| , |70| , 
for a general compact group it is constructed to 
all orders of PtQ. 

At present, chiral gauge theories obtained this 
way have a certain ambiguity which I would like 
to discuss briefly. A normalization factor (ab- 
solute value and phase) in Zi,, where is the 
partition function of the theory in the topological 
sector with charge v, is left undetermined. This 



is an unsatisfactory situation since, among oth- 
ers, it leaves the magnitude of fermion violating 
processes undetermined. 

In a recent paper||7^ H. Suzuki considered 
Weyl fermions in a real representation of the 
gauge group. He demonstrated that this Weyl 
fermion is equivalent to a Mayorana {L — R sym- 
metric) fermion on the lattice, like in the con- 
tinuum. Based on this observation Suzuki sug- 
gested a 'natural' relative normalization of differ- 
ent topological sectors for Weyl fermions in com- 
plex representations. 

The normalization of different topological sectors 
is related to CP symmetry in chiral gauge theories - 
a problem which keeps bugging me since 1998 ^ The 
present formulation is not explicitly CP symmetric. 
It might turn out that the renormalized final phys- 
ical predictions, which include the relative normal- 
izations between different topological sectors, respect 
CP. If this is the case, the existence of an explicit 
CP invariant lattice formulation is not excluded and 
would be interesting to find. If CP symmetry can not 
be restored (CP anomaly), that would be even more 
interesting. 

Under CP, ip -> IpW -W-'^tp ,U 

U'^^ , the vector gauge theory is invariant if 
WD{U'^^fW-^ = D(U), where C/^^ is the CP- 
transformed gauge configuration and in my conven- 
tion W = 72-^ I assume that D satisfies this con- 
dition. However, unlike in the continuum, the left 
handed part of the fermion action is not invariant: 

pp 1 

i^^D^PL^ iP^DiPl- -^D^r^Dip. (19) 

One might try to solve this problem by observing that 
the choice for the projection operators in eq. ([f7|) is 
not unique. One can find a continuous set of projec- 
tors constructed in terms of rg"-* = 75(1 — s_D) /A/"^"' 

and r^^' = (l-(l-s)D)75/A/''"' for the fermions and 
antifermions, respectively, where s is an arbitrary real 
parameter and M'^"^ is for the correct normalization. 
These projectors assure the L + R decomposition of 
eq. ( |l8| ) and give u for the difference between the di- 
mension of the left handed fermion and antifermion 

^Martin Liisclier raised my attention to the fact that the 
projection operators Pji/l in eq. ( |l7[ ) do not foilow the 
continuum transformation laws under CP. 

Under CP the spatial coordinates are reflected which we 
keep implicit in the following equations. 
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space^for any value of s. On the other hand, the CP 
violating term in eq. ( |l9| ) will be proportional to 1 — 2s 
and so for s = 1/2 the left handed action becomes CP 
invariant. Unfortunately, just at this value of s the 
projectors become singular (non-local). It seems to be 
difficult to define a CP invariant left handed action. 

Return back to the original definition in eq. jl^ ) 
and consider the fermionic expectation value on the 
gauge field configuration U 

{0)f = K,. J D^lD^^ eM-^LDiU)^L) O , (20) 

where O is some expression of the left handed fields 
and Kv is the unknown normalization factor of the 
topological sector u. Introducing a basis for the 
fermions and antifermions 

PlVj = Vj , ijjLix) = ^ CjVj{x) , (21) 

3 



eq. (pol) can be written as 



{0)f = K^YI J dcj J dck 



(22) 



exp(- 



k.j 



where Mkj = {u!k,Dvj). Going over to a different 
basis, the measure changes by a gauge field dependent 
phase. The difficult part of defining a chiral gauge 
theory is to fix this phase (or the basis) so that the 
final theory is local and gauge invariant ||6^-[7l[. 

The absolute value |(C)f| is, however, independent 
of the basis chosen. One might construct the {vj} 
and {iVk} basis vectors in eq. ( pl| ) in terms of the 
eigenvectors of D, for example. 

Consider a gauge field configuration U (U^^) with 
topological charge v = 1 (u = —1). D(U) has two left 
handed fermions with A = and 2, further a A = 2 
left handed antifermion. D{U'~^^) has a left handed 
antifermion with A = and no left handed modes 
with A = 2. The complex eigenvalues on U and U'^^ 
are the same. 

Let us chose a simple operator O — cq = (uq,^, 
where uq is the zero mode D{U)uq = 0. Eq. (p2) 
gives on the gauge field U 



\{0)f\ = \{co)f\^\K+,\2\{\X\ 



(23) 



where every complex pair of eigenvalues A, A* enters 
the product once. The factor 2 comes from the A = 2 
left handed mode. 

Compare eq. ( |2^ ) with 1(0) f| on the configuration 
U'^^ , where O = O'^^ = (%p,so) and so is the right 
handed zero mode on U'^^ . One obtains 



(24) 



In eq. (^|) the factor of 2 is missing: on U'^^ with 
V = —1 there is no left handed mode with A = 2. 

Suzuki, using a consideration unrelated to CP sym- 
metry, suggested a 'natural' relative normalization 
= 2""''^ This factor would restore CP sym- 

metry in the absolute value of the expectation values 
above. The question remains, however, whether this 
choice is consistent with cluster decomposition in the 
case where an instanton and anti-instanton are sepa- 
rated in a 1/ = configuration. 

4. Selected topics 

4.1. The topological susceptibility 

The topological susceptibility 



X 



d'^xq{x)q{Q)) 



V 



(25) 



where q{x) is the topological charge density and 
V = J d^xq{x). This quantity is related to one of 
the fanciest features of QCD, like the U{1) prob- 
lem, the axial anomaly and their interplay with 
topology. 

The topological susceptibility x"^ calculated in 
quenched QCD (i.e. over configurations gen- 
erated with the Yang-Mills action) enters the 
Witten-Veneziano relation |73| 



q _ ■I-!t"^v' 



2Nf ' 



(26) 



in the chiral limit. In full QCD, configurations 
with non-zero topological charge v are suppressed 
in the chiral limit by the fermion determinant. 
The leading correction to the chiral limit reads 



^This is the source of fermion number violation. 



X = 



Nf 



2Nf 



(27) 
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where — S = (uu) = (dd) ... is the quark conden- 
sate (in our normahzation f-^ w 93MeV in Na- 
ture). For a recent brief review on the history 
around these relations, see |7^, while for a recent 
summary on numerical simulations of I refer 
to [7^ . Here I would like to mention some recent 
developments only. 

On the theoretical side the GW formulation 
(which has exact SU{Nf) x SU{Nf) symmetry 
with the proper U{1) anomaly and the index the- 
orem) allows a clean derivation of eqs. (|6|j27|) 
on the lattice. Under the standard dynamical 
assumptions (chiral symmetry is spontaneously 
broken and the rj' remains heavy) eq. ( p7| ) was 
derived in|75|. In a recent paper the Witten- 
Veneziano formula was also reproduced on the 
lattice It is interesting to notice the analo- 
gies between the form of eqs. (^,^) and also in 
the steps of derivation. Although the 77' mass is 
due to a hard symmetry breaking (anomaly), it 
is switched off by the Nc — > 00 limit similarly as 
the soft breaking rUg is switched off in the chiral 
limit. 



0.08 



0.06 



0.04 



0.02 



0.00 



2 4 6 

Figure 5. Nf — 2 data: Wilson fcrmionsjs^ (open 
diamond), clover fermions|^(opcn circle), 
staggered fermions with smeared link opera- 
tor ||8^] (filled symbols), clover fermions with RG 
improved action[p8[ (burst). 



On the numerical side, the quenched topologi- 
cal susceptibility seems to settle down at a value 
X'^ = (0^05 MeV)'' with an estimated 10% error 
on X^ElJzll- ^ specialty of the topological sus- 
ceptibility is that its cut-off effects depend not 
only on the action but also on the details of the 
calculation (type and number of cooling steps) if 
non-chiral symmetric actions with a 'field theo- 
retic' definition of 1/ are used. The value above 
compares reasonably with eq. (|2^) if A^c = 3 num- 
bers are inserted on the r.h.s.. Nc = 3 seems to be 
close to the Nc ^ 00 limit which is also supported 
by a recent numerical study of the string tension, 
glue ball masses and x"^ in ^[/(TV), 2 < < SQ. 
The numerical analysis becomes cleaner if chiral 
symmetric actions are simulated. Only prelimi- 
nary results are available [ pOp^ . 

The topological susceptibility of QCD is, pre- 
sumably the best place to see, how sea quarks 
turn Yang-Mills gauge theory in full QCD. As 
eq. ( ^7| ) shows, the effect should be striking: 
the susceptibility x goes to zero in the chiral 
limit. The early results 1 81- 83 1 appear controver- 
sial. The authors in ref.[g5|86f argue that the cut- 
off effects are responsible for the confusing results, 
in particular, chiral symmetry violation for Wil- 
son type fermions and flavor symmetry violation 
for staggered fermions at large lattice spac- 
ings. The situation is expected to become better 
at smaller a, with improved actions, or improved 
topological charge operators. This expectation is 
confirmed by the Nf = 2 CP-PACS|8| and by 
the results in Fig. ||. 

Fig. ||is a compilation of staggered, Wilson 
and clover data with Nf = 2 at a « O.lfm, where, 
with the exception of the SESAM data|^, either 
the action, or the operator is improved. The ef- 
fect of dynamical fermions is seen clearly and the 
data are consistent with x ^ as ttItt — > 0. On 
the other hand, remaining discretization errors 
might be significant and the behavior of x with 
respect to (mTr^o)^ might look differently in the 
continuum limit ||87|. 

4.2. On quenched spectroscopy 

Quenched QCD is not a healthy QFT, but be- 
lieved to be universal in the continuum limit. The 
data are not very convincing, but at least do not 



10 



contradict this expectation. 



should be open minded. 



In my talk I referred to Sinya Aoki's contribution 
in Bangalore where the strong disagreement be- 
tween the staggered and Wilson Edinburgh plots was 
discussed. The situation concerning the hyperfine 
splittings confused me also. By now I understand 
these points better. 

'The staggered and Wilson fermion Edinburgh plots 
strongly disagree'. I believe, the reason of the dis- 
agreement in Fig. 4 inj89| are the very large cut-off 
effects in the unimproved staggered data which made 
a continuum extrapolation unreliable. This is nicely 
demonstrated in Fig. 19 of ref. [Q which shows that 
by reducing the cut-off effects in the staggered simula- 
tions (by simulating at small a, like injo^, or by using 
improved fermions) the staggered Edinburgh plot is 
pushed towards to that obtained with Wilson type of 
fermions. 

'Are the quenched hyperfine splittings too small, or too 
large?' They were found too small by many groups 
and most convincingly by the CP-PACS CoUabora- 
tion[Q . On the other hand, UKQCD finds the hyper- 
fine splittings always far too large . I think the rea- 
son is that UKQCD fixed the scale from ro ~ 0.5fm, 
while the standard choice is fixing mp. Fixing rup 
gives ro ~ 0.55fm (see, for example [^) and this 10% 
effect explains the paradox with the hyperfine split- 
tings. The disagreement of quenched spectroscopy 
with Nature can be shifted to different quantities by 
changing the method of fixing the scale. 

4.3. Strange cut-off effects 

Numerical data in asymptotically free theories 
are standardly extrapolated to the continuum us- 
ing Symanzik's results |2^] valid in every order of 
perturbation theory. Accordingly, in a bosonic 
theory the leading cut-off effect is quadratic O(a^) 
modified by logarithms. Results in the 0(3) 
non-linear sigma model show that either non- 
perturbative effects invalidate this expectation, 
or the onset of leading behavior happens at very 
large correlation lengths only|Q. In Fig. ^ the 
data obtained with two different actions on the 
step scaling function follow a linear 0{a) type 
of behavior - at least up to correlation lengths 
350. The cut-off dependence of the zero mo- 
mentum 4-point function is similar pM|j95[|. One 



1.31 
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0.05 
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Figure 6. The step scaling function 'S{2,uo,a/L) 
at uq — 1.0595 for two different actions. 
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